Isobaric multiplet yrast energies and isospin non-conserving forces 
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The isovector and isotensor energy differences between yrast states of isobaric multiplets in the 
lower half of the pf region are quantitatively reproduced in a shell model context. The isospin non- 
conserving nuclear interactions are found to be at least as important as the Coulomb potential. Their 
isovector and isotensor channels are dominated by J = 2 and J = pairing terms, respectively. The 
results are sensitive to the radii of the states, whose evolution along the yrast band can be accurately 
followed. 

PACS numbers: 21.10.Sf, 21.60.Cs, 23.20.Lv, 27.40.+Z, 29.30.-h 



The electrostatic energy of a sphere of radius R and 
charge Ze is easily calculated to be Ec = 3e 2 Z 2 /5R. It is 
under this guise that the Coulomb field enters the Bethe 
Weizacker mass formula, and becomes a basic quantity 
in nuclear structure. Direct evidence of entirely Coulomb 
effects has long been available from displacement ener- 
gies between mirror ground states (MDE), and in the 
last decade from differences in excitation energies of yrast 
bands in mirror nuclei (MED) @, §, | [§ §|(J. 

The MDE energies range from few to tens of MeV. 
They are given mainly by Ec, but precise calculations 
were found to be unexpectedly inaccurate — the Nolen- 
Schiffer anomaly || — and revealed the necessity to in- 
troduce charge symmetry breaking (CSB) nuclear poten- 
tials (j7| and references therein). The anomaly is now 
under control to within shell effects (§1, [|, which we de- 
fine as deviations from a Bethe Weiszacker-type formula 
involving only number of particles (A), isospin (T), and 
its third component (T z ). 

The MED are defined by (Z > and Z < are the largest 
and smallest Z in the multiplet, and Ej are the yrast 
excitation energies) 



shell): 



/#> = VbI - V£, = VIZ + VZ - 2VSr- (3) 
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nil in Vb , while for Vc we have /3q1 = (3^1 

The MED are entirely of isovector origin and the first 
exact shell model calculations in the full pf shell indi- 
cated that Vq° , i.e., calculated in the harmonic oscillator 
(ho) basis fails to give a satisfactory description || . The 
way out proposed in this reference consisted in replacing 
the harmonic oscillator matrix elements Vil, by em- 

Ot 7 / 2 J 

pirical ones derived from the A = 42 spectrum which are 
very different. Therefore, it was hard to attribute the 
replacement to a renormalization of Vc, expected to be 
small. But it was equally hard to think in terms of CSB 
precisely because the effect was so large. Nevertheless the 
ansatz (or variants of it) worked quite well, and subse- 
quent calculations incorporated it ||, |j leading eventu- 
ally to (almost) full quantitative agreement |L2| for the 
MED in A = 47, 49, 50 and 51. When the isotensor 
TED data came in, it became clear that both charge in- 
dependence breaking Q and CSB had to be invoked. 



MED j = Ej(Z + T) — Ej(Z -T), Z = 



(1) 



The observed MED are very small (of the order of 10- 
100 keV), and entirely due to shell effects. Recently, the 
experimental information on yrast bands has been ex- 
tended to isospin triplets |l^, [ll|] , thus determining new 
quantities, the TED given by 

TED j = Ej(Z + 1) + Ej(Z - 1) - 2Ej(Z). (2) 

Both measurements are needed to achieve a clear under- 
standing of the interplay between the Coulomb potential 
Vc, and Vb, the isospin breaking nuclear interaction. 

To analyze them, we start by writing the isovector, 
pi 1 ^, and isotensor, [3 r 2 \ contributions to Vb, as linear 
combinations of two body matrix elements in neutron- 
proton (f7v) formalism (r = r% r2 ^3 r^, where fj is a sub- 



TABLE I: Coulomb (Vb), isovector (MED-Vb 

isotensor (TED-Vb 
calculated in the oscillator basis (ho) 



/3f 2 ' ) energies (keV) in A 
7/2 



and 
42. Vb 



J = J = 2 J = 4 J = 6 



V c - 
£M 42 Ti 



v Cf. 



Sj [ 42 Ti + 42 Ca 



7/2 

Ca] - Vc 
- 2 42 Scl- 



Vc 



81.60 
5.38 
116.76 



24.60 
92.55 
80.76 



6.40 -11.40 
4.57 -47.95 
2.83 -42.15 



Ironically, this result is obvious from the, long known, 
A = 42 spectra [B. Assuming that the observed states 
are essentially fy 2 configurations on top of the 40 Ca 
core, these spectra define an interaction in the sub- 
shell. Therefore, by setting Vc = ^ / cf 7/2 7 the nuclear 
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isovector and isotensor contributions can be extracted. 
They are shown in Table |, where their centroids - 
J2j( 2J + 1 ) v / r/2 / T,j( 2J + 1 ) have been subtracted for 
clarity. The assumption of fy 2 dominance is not war- 
ranted, as — at least — the J = and 2 states are known 
to mix with core excitations. Therefore a safer procedure 
consists in replacing the lowest observed states by the 
fy 2 centroids estimated from spectroscopic factors [jli] . 
However, when this is done, no significant change obtains 
in Table |[ whose indications must therefore be taken very 
seriously. 

A renormalized Vc adapted to the fy 2 space will re- 
main a purely itit force, and therefore the same for the 
isovector and isotensor channels. Furthermore, it is ex- 
pected to be reasonably close to the bare Vc in the first 
line of Table |. Therefore, upon subtracting this bare Vc 
from the observed data — second and third lines — we ex- 
pect the same, reasonably small numbers. It is obvious 
that the corresponding numbers are neither equal nor 
small. The unavoidable conclusion is that the A = 42 
data indicate that the role of isospin non conserving nu- 
clear forces is at least as important as that of the Coulomb 
potential in the observed MED and TED. 

For the full description of these quantities in A = 46-51 
we rely on exact, isospin conserving shell model calcula- 
tions pl| with single particle spectrum from 41 Ca and the 
KB3G interaction. Very little changes are observed if the 
other standard interactions are used (KB3, FPD6, all de- 
fined in [^6|). The energy differences are obtained in first 
order perturbation theory pl[ , as the sum of expectation 
values, in which we separate the monopole and multipole 
components of the Coulomb field Vc = Vc m + Vcm fol- 
lowing Refs. §, ||: 

MED j - A M {Vc m )j + A M {VcM)j + A M {VB)j,(4) 
TED,/ = A t (V C m)j + A t (V b )j. (5) 

The monopole Vcm contains all terms quadratic in scalar 
products of Fermion operators af ■ aj . The non diagonal 
contributions (i =^ j) lead to isospin mixing that demands 
second order perturbation theory. They will be consid- 
ered here only through their influence on the radial wave- 
functions, i.e., the Thomas Ehrman shift that depresses 
the single particle p 3 / 2 state in 41 Sc by 200 keV below its 
analogue in 41 Ca. The diagonal part (i — j) involves only 
proton number operators. It contains Ec plus a single 
particle splitting induced by Vc on the orbits of principal 
quantum number p above harmonic oscillator (ho) closed 
shell Z r 
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-0.001 Z, 



13/12 



[21(1 + 1) - p(p + 3/2)] 



AV3(p + 3 /2) 



MeV. (6) 



The effect of Ec is proportional to the difference of (in- 
verse) radii between a J-yrast and the ground state [0. 
The total radii depend on those of the individual orbits, 



and therefore — to good approximation — on the average 
neutron plus proton occupancies for each orbit, which we 
denote by (mk) j/2, with mk = Zk + (number of neu- 
trons plus number of protons in orbit k). We take aver- 
ages relying on the near equality of proton radii in both 
members of a mirror pair @. As it is reasonable to as- 
sume that orbital radii depend only on I, and the P1/2 oc- 
cupancy is always negligible, the whole radial effect will 
be taken to depend on the p 3 / 2 occupancy. Note that the 
single particle contribution from Eq. (Q) is proportional 
to the difference of proton and neutron occupancies. It 
is important in A = 41, but typically ten times smaller 
than the radial effect in A =47-51, so we neglect it and 
end up with A M (Vcm)j = a m {m P3/2 ) j /2. The value of 
a m can be estimated by adding to the observed shift the 
single particle splitting @ that depresses the I = 3 orbits 
with respect to the / = 1 ones (by 125 keV at A = 40, 
Z cs = 20). Then, a m w .200 + .125 = 0.325 MeV. 

In the isotensor case the rn P3/2 contributions cancel 
out. 

The multipole contribution A(Vcm)j is calculated us- 
ing oscillator Coulomb matrix elements in the pf shell. 
The only direct information on Vb comes from Table @. 
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FIG. 1: Example of MED renormalization in A = 50 



To make use of it we must explore the possibility of spec- 
ifying an interaction acting in the full pf shell, solely in 
terms of j V/2 matrix elements. The idea turns out to be 
quite viable using the multiplicative prescription J23[ 



A(V$; t )j = bA(V£ )j, 



(7) 



as illustrated for the Coulomb potential in Fig. [l| 

The same form efficiently relates the schematic pair- 
ing or quadrupole pairing forces in the pf shell to the 

yj=0or j=2 ma ^ r j x elements. To minimize the num- 
17/2 

ber of parameters, for Vb we retain only the leading 
term suggested by Table |, and set 8^ — 0i , 

B^ 2 ) — do V/~°, where V/ is the matrix element with 

^pi ^ t 7 / 2 ' I7/2 

unit value. Collecting all the pieces we have 
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FIG. 2: Experimental §, |, g, | and calculated MED for the pairs 47 V- 47 Cr, 49 Cr- 49 Mn, 50 Cr- 50 Fe, and 51 Mn- 51 Fe. 



MED j = -a m {m pa/2 )j + A M {v£° f + [3, V f ' T = 2 )j, (8) 

TED j = A T (v£° pt + f3 2 V f J 7 =°)j. (9) 

In Figs. H for the MED, Vc m , Vcm, and Vbm stand 
respectively for the first, second and third terms in 
Eq. (||). The parameters are taken to be round numbers, 
a m = 300 keV and ft = ft = 100 keV. 
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FIG. 3: TED for A = 46 and 50 

The reduction of Vb, for MED and TED, to a single 
matrix element is an oversimplification, but the results 
are so satisfactory that the need of extra terms is not felt. 



The only parameter-free alternative is to take matrix el- 
ements with the weights in Table |j. However, this choice 
is arbitrary because — from the discussion around Eq. (0) 
and Fig. |l| — we expect a case by case (even matrix ele- 
ment by matrix element) renormalization. Nonetheless, 
though the agreement with experiment becomes less im- 
pressive, it remains acceptable. The conclusion is that 
the leading term in J = 2 for MED is indeed dominant, 
and that in J = for TED very dominant. 

The Vcmj Vcm and Vbm contributions, shown sepa- 
rately in Fig. H for A = 47 and 49 are quite far from the 
observed pattern, which is accurately reproduced only 
after these disparate terms are added. For A = 50 and 
51 we have replaced the Vbm part by a variant of the 
full sum in which ft is halved. For A — 50 the changes 
are insignificant, but there is a definite improvement in 
A = 51 (remember again Eq. (fjj) and footnote before it). 

It is especially worth noting in Fig. || that the strong 
signature effect in the A = 49 band is erased in the MED 
by the out-of-phase Vc m and Vcm, while the signature 
staggering is enhanced in A = 51. 

The experimental TED patterns in Fig. |^ are quite 
nicely reproduced by the minimal = ft V/^j® choice. 
As mentioned, the inclusion of the J ^ terms (third 
line of Table |) makes little difference, and — interestingly 
enough — simply ignoring Vb and doubling Vcm (or the 
other way round) makes practically no difference. Which 
confirms the overwhelming dominance of J = pairing. 
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It can be hoped that a rigorous treatment calling upon 
state of the art CSB potentials will confirm the role 
of the J = 2 pairing term for the isovector MED. The 
TED behaviour seems far simpler and our results are 
consistent with the findings in for (3^ borne out 
in |ll[] for A = 46. Therefore, here we may bet on — rather 
than hope for — confirmation by the charge independence 
breaking potentials [O . 

The isovector channel raises a difficulty for A = 46. 
In 0, |l8f it was found that w using the same 

(2) 

functional form as for p). with strong J = pairing, 
which does not square with our results. But in this case 
our results do not square with experiment either. The 
scheme that has been successful in A = 47, 49, 50 and 51 
fails in A = 46: we are simply unable to do any better 
than in |I| Fig. 3a)]. 



that direct evidence for charge symmetry breaking has 
been confined, so far, to the very light systems (basically 
A = 2 and 3) 0. The mechanism plays an important 
part in resolving the Nolen Schiffer anomaly in the MDE, 
but the effects of Vc remain much stronger^. For the 
MED and TED, Vc is at most as strong as Vb, for which 
we have shown that substantial quantitative information 
can be extracted from the data. To boot, the MED also 
provide a view of the evolution of yrast radii. 

This work owes much to a stay of AZ at the UAM, 
made possible by a scholarship of the BBVA founda- 
tion. AP is supported by grant BFM2000-30 from MCyT 
(Spain). 




FIG. 4: Yrast energy differences in A = 46 



The trouble is no doubt due to the poor spectroscopy 
provided by full pf shell diagonalizations for A < 46, 
at least when compared with the very high quality de- 
scriptions for the rest of the f?/ 2 nuclei (i.e., A < 56). 
Fig. |] illustrates the point: the calculated yrast ener- 
getics is wrong for the lowest states and, for the others, 
far less precise than the corresponding patterns in the 
heavier nuclei (see for example jl9) for A — 47 and 49). 
This problem extends to transition rates and static mo- 
ments. It was first noted and abundantly discussed in 
Rcf. p] ] but its quantitative explanation remains a chal- 
lenge. This unsatisfactory situation provides nonetheless 
a helpful clue: the TED may be unsensitive to details, 
but the MED demand accurate wave functions and could 
be taken as tests of their quality. 

Within the A = 46 proviso, our results make obvious 
something that may seem at first surprising: isospin non 
conserving potentials play a role that is at least as impor- 
tant as Vc in explaining the MDE (and TDE, as found 
previously in [fy]). In this respect, it is worth noting 
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